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Hadronic B decays 

M. Beneke a * 

a Institut fiir Theoretische Physik E, RWTH Aachen, Sommerfeldstr. 28, D - 52074 Aachen, Germany 

I briefly summarize the factorization approach to hadronic B decays emphasizing theoretical results that have 
become available recently. The discussion of its application to data is abridged, and only the determination of 
7 = (71 ± 5)° from time-dependent CP asymmetries is included in some detail. 
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1. INTRODUCTION 

Many observables at the B factories are con- 
nected with branching fractions, CP asymmetries 
and polarization of exclusive, hadronic B decays. 
They provide access to the flavour and spin struc- 
ture of the weak interaction, but a straightfor- 
ward interpretation is usually obscured by the 
strong interaction. In technical terms, the diffi- 
cult (long-distance) part of the strong interaction 
resides in the matrix elements (f\Oi\B), where 
Oi is an operator in the effective weak interac- 
tion Lagrangian. 

Systematic approaches to hadronic B decays 
are based on expansions in small parameters. The 
two available options exploit approximate flavour 
symmetries (expansion parameter m q /A, m q a 
light quark mass), or the large energy transfer 
in B decays (expansion parameter A/rrib), result- 
ing in two frameworks - "SU(3)" and "Factor- 
ization" - that could hardly be different method- 
ically and technically. In practice, both frame- 
works are implemented only at the leading order, 
and additional assumptions are usually necessary 
(neglecting "small" amplitudes; estimating K/rrib 
corrections). Despite this restriction, there has 
been much progress by applying and working out 
these theories over the past few years. In the fol- 
lowing I focus on the factorization approach. Fur- 
thermore, / will be assumed to be a charmless, 
two-body, meson final state; the mesons are as- 
sumed to be pseudoscalar or vector mesons from 
the ground state nonet. 
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2. THEORY OF HADRONIC DECAYS 
(FACTORIZATION) 

The starting point is the investigation of Feyn- 
man diagrams with external collinear lines (en- 
ergetic, massless lines with momenta nearly par- 
allel to one of the two final state mesons, Mi or 
M2), one nearly on-shell heavy-quark line, and 
soft lines (representing the light degrees of free- 
dom in the B meson). The simultaneous rele- 
vance of collinear and soft configurations implies 
three relevant scales: m&, y/rribA, and A. In the 
heavy-quark limit the first two are perturbative, 
and only the third is long-distance. Factorization 
amounts to showing that the long-distance contri- 
butions to the matrix elements ([cl] [c2]|Oj| [s]) are 
actually contained in the simpler matrix elements 
([cl]|(g&)(0)|[*]> (form factors), ([ci]|g(aj)g(0)|0>, 
and (0|g(x)&(0)|[s]) (light-cone distribution am- 
plitudes). It is then assumed that if this holds 
perturbatively to all orders for all quark-gluon 
matrix elements, then it does for the hadronic 
matrix elements. 

Factorization in a similar form was first applied 
to B decays in [1] as a phenomenological approx- 
imation akin to the vacuum saturation approx- 
imation for the four-quark operator matrix ele- 
ments relevant to BB mixing. Intuitively, fac- 
torization might work, because the partons that 
eventually form the meson M% that does not pick 
up the spectator quark escape the B remnant as 
an energetic, low-mass, colour-singlet system, and 
hadronize far away and therefore independently 
from the remnant. This qualitative argument was 
given in [2] for the decay Bd — * D + tt~. In [3] it 
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was shown to hold for charmless decays, where 
the disruption of the B meson is much more vio- 
lent, and a calculational framework was provided, 
in which the phenomenological factorization ap- 
proach was contained as a leading-order approx- 
imation. At the same time, the next-to-leading 
order corrections were computed. 

The new factorization formula included a new 
mechanism, spectator-scattering, where a hard- 
collincar interaction with the soft remnant takes 
place. Thanks to the development of soft- 
collinear effective theory (SCET), this mechanism 
is now much better understood. In the following 
I sketch the rederivation of the factorization for- 
mula in SCET [4-6]. 

Integrating out fluctuations on the scale m& at 
leading power in the A/mj expansion amounts to 
an analysis of the structure of hard subgraphs 
with external hard-collincar, collinear and soft 
lines. Those identified as leading are then calcu- 
lated perturbatively in a s (nib). Formulated as an 
operator matching equation from QCD to SCETi, 
the result of this analysis reads 

+ Cj 1 (t, S )*[t4Asn+)h v ]y (1) 

Remarks: (a) The short-distance coefficient C\ 
incorporates corrections to naive factorization. 
The term in the second line describes spectator- 
scattering with its own short-distance coefficient 
Cj l . (b) The second line is a leading contribu- 
tion despite the fact that the corresponding op- 
erator is suppressed in dimensional and SCETi 
power counting. This follows by extension of the 
power-counting analysis of [5] . (c) The meson M2 
factorizes already below the scale nib [6], since 
SCETi does not contain interactions between the 
X*- -* fields and the collinear-1 and soft fields. It 
follows that at leading power in the heavy-quark 
expansion, the strong interaction phases originate 
from the short-distance coefficients C\' 11 at the 
hard scale, (d) The result above must be modi- 
fied to account for a non-factorizing effect when 
the final state contains an rf'^ meson. This ef- 
fect is explained in [7], but appears to have been 
missed in the SCET rederivation of the factoriza- 
tion theorem for mesons with flavour-singlet com- 



ponents [8]. Taking the hadronic matrix element 
of (1) gives 

(M 1 M 2 |O i | J B) = $m 2 N* (t\u)F bm ^0) 

+ C u (t,u)*E bm ^t,0)) 7 (2) 

where I have reintroduced the full QCD form 
factor F SMl (0) resulting in a slight modifica- 
tion of the short-distance coefficients. E BMl (r, 0) 
denotes a new, unknown, non-local form factor, 
which depends on the convolution variable r. 

The different implementations of factorization 
can be distinguished broadly by their treatment 
of the different factors in (2). In the PQCD 
approach [9] the form factors F BMl (0) and 
E BMl (t, 0) are assumed to be short-distance dom- 
inated, and claimed to be calculable in a gener- 
alized factorization framework (fc^-factorization). 
All four quantities, T 1 , C 11 , F BMl (0), (t, 0) 
have been calculated at leading order. Recently, 
some next-to-leading order (NLO) corrections to 
T 1 have been included. In the QCD factoriza- 
tion approach [3] it is assumed that the stan- 
dard heavy-to-light form factors receive a lead- 
ing soft contribution, and are therefore not cal- 
culable. However, S BMl (r, 0) is dominated by 
pcrturbative hard-collincar interactions, and fac- 
torizes further into light-cone distribution ampli- 
tudes (see below). In the BBNS implementation 
of QCD factorization, F BMl (0) is a phenomeno- 
logical input (usually from QCD sum rules). The 
other three quantities, T l , C u , ~ BMl (r, 0) have 
been calculated at the next-to-leading order. In 
the BPRS implementation [6] the use of pertur- 
bation theory at the hard-collinear scale y/nibA is 
avoided, and both form factors are fit to hadronic 
B decay data. This approach is restricted to 
leading-order in the short-distance coefficients, 
since only then does the unknown form factor 
E BMl (r, 0) enter the equations through a single 
moment. There is another difference between 
BBNS and BPRS, who claim that (1) is not valid 
for diagrams with internal charm quark loops. 
(This should be distinguished from [10], which 
speculates about large power corrections from 
charm loops or annihilation.) I believe that the 
theoretical arguments leading to this conclusion 
are wrong [11]. For phenomenology, the impor- 
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tant consequence from treating charm loop di- 
agrams as non-perturbative is that the penguin 
amplitudes must be determined from data, such 
that no CP asymmetry can be predicted from 
theory alone. Since the tree amplitudes are also 
determined from data (namely, through the two 
form factors; the phase of C/T is automatically 
zero in a leading-order treatment), the BPRS ap- 
proach has much more in common with amplitude 
fits to data than with QCD/SCET calculations. 

The QCD factorization argument is completed 
by noting that the non-local SCETi form factor 
S BMl (t, 0) factorizes into light-cone distribution 
amplitudes, when the hard-collinear scale ^/nibA 
is integrated out [5]. Inserting 

S BMl (r,0) = J(t;w,v) * $ B M * $ M » (3) 

into (2) results in the original QCD factoriza- 
tion formula with the additional insight that the 
spectator-scattering kernel T u = H u * J factor- 
izes into a hard and hard-collinear kernel. The 
development of SCET was crucial to identify the 
operators and precise matching prescriptions that 
make the calculation of higher-order corrections 
to spectator-scattering feasible. 

3. HIGHER-ORDER CALCULATIONS 

On the calculational side one of the main ef- 
forts over the past few years has been the cal- 
culation of one-loop corrections to spectator- 
scattering, which formally represents a next-to- 
next-to-leading contribution in the QCD factor- 
ization approach. This programme is now com- 
plete. The hard-collinear correction to J has been 
calculated in [12-14]; the hard correction to H u 
in [15,16] for the tree amplitudes and in [17] for 
the QCD penguin and electroweak penguin am- 
plitudes. (An earlier calculation of the QCD pen- 
guin contribution [18] disagrees with [17].) The 
main results are summarized as follows: (a) The 
convolution integrals are convergent, which estab- 
lishes factorization of spectator-scattering at the 
one-loop order, (b) Perturbation theory works for 
spectator-scattering, including perturbation the- 
ory at the hard-collinear scale, (c) The correc- 
tion enhances the colour-suppressed tree ampli- 
tude, and reduces the colour-allowed one. This 



improves the description of the tree-dominated 
decays to pions and p mesons, (d) The correction 
to the colour-allowed QCD penguin amplitude 
is negligible. Thus there is no essential change 
in the predictions of branching fractions and CP 
asymmetries of penguin-dominated decays. 

The evaluation of the colour-suppressed tree 
amplitude gives [17] 

a 2 (7T7r) = 0.18- [0.15 + 0.08i] N LO 

sp 1 {[0.12] LO + [0.05 + 0.05z]nlo 

+ [0.07] tw3 }. (4) 
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Here r sp = (9f Ml fB)/(m b F BM ^0)X B ) defines 
a combination of hadronic parameters that nor- 
malizes the spectator-scattering effect. Eq. (4) 
shows the importance of computing quantum 
corrections: the naive factorization value 0.18 
is nearly canceled by the 1-loop vertex correc- 
tion calculated in [3]. It now appears that the 
colour-suppressed tree amplitude is generated by 
spectator-scattering. It is not excluded that r sp 
is a factor of two larger than 0.485, in which case 
ai becomes rather large. My interpretation of the 
pattern of the irir, np and pp branching fractions 
is that spectator-scattering is important [19]. On 
the other hand, the large direct CP asymmetry 
in B d — > 7r + 7r~ cannot be explained by known 
radiative corrections, and remains a problem. 

Next-to-leading order corrections have recently 
been implemented in the PQCD approach for the 
first time [20]. More precisely, the 1-loop ker- 
nel T l from the QCD factorization approach is 
used as a short-distance coefficient for the sub- 
sequent tree-level PQCD calculation. The nu- 
merical impact is again strongest on the colour- 
suppressed tree amplitude, C. But while this cor- 
rection (— [0.15+0. 08i]nlo m (4)) results in a near 
cancellation of the naive factorization term in the 
QCD factorization approach, it provides an en- 
hancement of C by a factor of several in [20] . This 
resolves the nK puzzle in the PQCD approach. 

I am rather sceptical about the possibility to 
perform accurate calculations in the PQCD ap- 
proach. A complete NLO calculation in the 
PQCD approach requires a calculation of all 
one-loop spectator-scattering diagrams (similar 
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to [15-17]) rather than the 1-loop BBNS kernels. 
The calculation of T 1 is done with on-shell exter- 
nal lines, but when the vertex diagram appears 
as a subdiagram in a larger diagram with hard- 
collinear exchanges, the external lines of the sub- 
diagram can be far off-shell. Hence T 1 is not the 
appropriate quantity to be used. The numerical 
differences between (4) and [20] despite the same 
input T 1 can be traced to the choice of scales. 
The one-loop correction to T 1 makes the result 
less sensitive to variations of the renormalization 
scale in the Wilson coefficients, but only for scales 
larger than about 1.5 GeV, below which pertur- 
bation theory breaks down. Factorization shows 
that the scale of the Wilson coefficients should be 
of order nib- However, in the PQCD approach 
the scales m,b and V m bA are not distinguished, 
and the Wilson coefficients are evaluated at very 
low scales (to 500 McV) , where perturbation the- 
ory is not reliable. An unphysical enhancement 
of the Wilson coefficients at small scales is also 
the origin of the large penguin and annihilation 
amplitude in the PQCD approach. Yet a varia- 
tion of the renormalization scale is not included 
in theoretical error estimates. 

4. POWER-SUPPRESSED EFFECTS 

Power corrections to the QCD penguin ampli- 
tudes are essential for a successful phenomenology 
within the factorization framework. The most im- 
portant A/nib effect is the scalar QCD penguin 
amplitude r x a 6 . Fortunately, the bulk contribu- 
tion to this amplitude appears to be calculable, 
although its factorization properties are not yet 
understood. This power correction is responsible 
for the differences between PP, PV and VV fi- 
nal states and the rj^K^ final states [7]. The 
calculated pattern is in very good agreement with 
experimental data. 

The second most important power correction 
is presumably weak annihilation. I emphasize 
"presumably" , since there is no unambiguous em- 
pirical evidence of any weak annihilation con- 
tribution in charmless decays, and only upper 
limits can be derived. The theoretical difficulty 
with power corrections is reflected in the different 
treatments of annihilation. In PQCD it is calcu- 



lable and large. In the BBNS implementation of 
factorization it is represented by a phenomenolog- 
ical parameter [21], not very large, but it makes 
the calculation of CP asymmetries uncertain. In 
the BPRS implementation it is neglected together 
with all power corrections. This is phcnomcnolog- 
ically viable, since the charm penguin amplitude 
is fit to data anyway. Some weak annihilation 
amplitudes have been calculated with light-cone 
QCD sum rules [22]; the result is compatible with 
the BBNS parameterization. 

It is not difficult to write down the power- 
suppressed operators in SCET [23]. The problem 
is that the factorization formula involves convo- 
lutions, which usually turn out to be divergent 
at the endpoints, making the result meaningless. 
The inadequacy of SCET in addressing this well- 
known problem in hard-exclusive scattering was 
pointed out in various forms in [5,24], but no so- 
lution was offered. In the recent paper [25] it is 
proposed that endpoint divergences can be elimi- 
nated by a new type of factorization ( "zero-bin" ) . 
This would be a breakthrough; however, I do not 
see how "zero-bin" factorization could possibly 
be correct, since it cuts off the endpoint con- 
tributions without defining the appropriate non- 
pcrturbative objects that would represent the 
endpoint region. Thus, the new factorization- 
scale dependence is not consistently canceled. 

To explain this I compare the treatment of a 
certain weak annihilation diagram in "zero-bin" 
factorization [26] with the BBNS parameteriza- 
tion [19,21]. In the first method, the divergent 
integral a s J Q dx 4>m 2 i x ) /x 2 is interpreted as 



-<// M2 (l).a s ln^ + F, 

M— 



(5) 



in the second as 

-<f>' Ma (l)(l + g A e^) -a s ln^ + F, (6) 
where 

<j)M 2 {x) +X(j>' M2 (l) 



n, I dx ■ 

'o 



(7) 



is a finite, subtracted integral. In [26] [i- is taken 
to be of order m^, thus the first term in (5) is of 
order a s and perturbative. The endpoint contri- 
bution is effectively set to zero, but the depen- 
dence on the arbitrary factorization scale yu_ is 
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not canceled. A candidate non-perturbative pa- 
rameter for the endpoint contribution could be 
(f>' M2 (l), but this object is not defined in SCET, 
so a field-theoretical definition of the method is 
missing. The second expression (6) looks similar, 
but now there is a large endpoint logarithm, and 
a s lnms/A is of order 1. The endpoint contribu- 
tion is considered to be non-perturbative, and is 
parameterized by the complex quantity qa& iva ■ 
It is again the absence of a field-theoretical defi- 
nition of this quantity that makes the BBNS pa- 
rameterization a phenomenological model. Ex- 
pression (6) is clearly a more conservative treat- 
ment of the problem than (5). 

It is evident that in the absence of a field- 
theoretical definition of the zero-bin subtraction 
method, the statement that "annihilation is real 
and calculable" is wishful thinking (I share the 
wish.); it also contradicts the QCD sum rule cal- 
culation [22]. My strong criticism (prompted by 
strong claims) is not to mean that the problem 
of endpoint factorization is not important. To 
the contrary, its solution is prerequisite to further 
progress in SCET. 

5. PHENOMENOLOGY (OMITTED) 

There is not enough space to discuss the factor- 
ization calculations of branching fractions and CP 
asymmetries and the comparison with data. I fo- 
cus on the calculation of the CP-violating S pa- 
rameters and the determination of 7 in the fol- 
lowing section. A very brief summary of the other 
topics discussed in the talk reads: 

• The global comparison of all B — ► PP, PV 
data with scenario S4 of [19] remains im- 
pressively good, including CP asymmetries, 
but there are persistent exceptions. The 
same is true for the PQCD [20,30] and 
BPRS [8,31] approaches. 

• An enhancement of the electroweak pen- 
guin amplitude to explain the nK system 
is no longer compelling. The difference be- 
tween the CP asymmetries in ir°K ± and 
ir^K seems to require an enhancement 
of the colour-suppressed tree amplitude, 
which cannot be explained by factorization. 



• There exist interesting effects [27,28] in 
B — ► VV decays, which motivate further 
polarization studies. See [29] for a compre- 
hensive analysis of these decays. 
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Figure 1. CKM phase 7 from Sf with / = 

(7T7T,7rp, [pp] L ). 

6. DETERMINATION OF 7 FROM S f 

The time-dependent CP asymmetries Sf in tree- 
dominated AD = 1 decays are particularly suited 
[19,21] to determine the CKM phase 7 (or a; I as- 
sume that is determined experimentally) in the 
framework of QCD factorization, since hadronic 
uncertainty enters only in the penguin correction; 
the dependence on strong phases is reduced, be- 
cause it arises only through cos<5; the sensitivity 
to 7 is maximal near 7 ~ 70°. 
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Figure 1 shows that for / = (tttt, np, \pp\i) 
and measurements S f = (-0.59 ± 0.09, 0.03 ± 
0.09, -0.06 ±0.18) (HFAG averages), one obtains 
(ignoring a second solution that does not lead to 
consistent results) 7 = (70±^,69 ± 7, 73 ± 8)°. 
The three determinations are in agreement with 
each other, resulting in the average 7 = (71 ±5)°. 
See [19,29] for details. 



7. CONCLUSION 

The subject of hadronic decays has been and still 
is a very fertile ground for developing new the- 
oretical concepts in heavy flavour physics. A 
lot has been learned about hadronic dynamics. 
Moreover, 7 is by now known rather well from 
charmless decays. There should be some way to 
include this information in the CKM fits. 
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